On 
On 



Generalised Functions and Distributional Curvature 

of Cosmic Strings 

C J S Clarke, J A Vickers and J P Wilson 
Department of Mathematics, University of Southampton, Southampton, S017 1BJ, UK 

Abstract A new method is presented for assigning distributional curvature, in an invariant manner, 
to a space-time of low differentiability, using the techniques of Colombeau's "new generalised functions" . 
The method is applied to show that curvature of a cone is equivalent to a delta function. The same is 
true under small enough perturbations. 

1. Introduction 

For some time it has been recognised that there is an important place in relativity for metrics whose 
curvature has to be regarded as a distribution; that is, the components of the Riemann tensor have to be 
interpreted not as ordinary functions, but as functionals (Gel'fand and Shilov, 1963) or as ideal limits of 
qq , functions. This has been applied successfully to such examples as surface-distributions of matter (Israel, 

£S) ' 1966; Clarke and Dray, 1987) and to gravitational radiation (Kahn and Penrose, 1971). In these cases, 

however, it is possible to formulate the field equations so as to avoid operations that are not well defined in 
ordinary distribution theory: specifically, multiplying distributions, or multiplying, say, a delta-function 
by a discontinuous function. The potential ambiguity of such operations was established by Schwarz 
(1954), who pointed out that it was impossible consistently to define an associative multiplication on 
distributions, together with an operation of differentiation, which coincided with the usual definition of 
these on continuous and C 1 functions, respectively, and for which there was a non-zero distribution S 
, satisfying x5(x) = 0. 

The curvature tensor is a non-linear function of the metric, so that if one is to avoid these illegitimate 
operations, then there are strong constraints on the sort of metric that can be considered. Specifically, 
Geroch and Traschen (1987, Theorem 1) showed that in order for the components of the Riemann tensor 
and its contractions to be well defined as distributions, their singular parts have to have support on a 
submanifold of dimension of at least three. This unfortunately excludes many space-times which one 
expects to have a distributional curvature, for example the cone-like space-time 

ds 2 = -dt 2 + dr 2 + A 2 r 2 d(/) 2 +dz 2 , \A\ < 1. (1) 
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Various attempts (Raju, 1982; Balasin and Nachbagauer, 1993) have been made to give prescriptions 
for multiplying distributions in order to overcome the problem posed by the nonlincarity of Einstein's 
equations, but they have tended to depend on regularization procedures whose invariance and general 
applicability were uncertain. There has been until recently no general theory within which it could be 
shown that the result was invariantly defined, independently of the particular regularization adopted. 

Our aim in this paper is to describe within the setting of relativity theory the formalism, due 
to Colombeau (1983, 1990), in which distributional curvatures can be rigorously and unambiguously 
defined. We will illustrate the method in the case of the conical metric (1), where one expects there to 
be a distributional energy momentum tensor of the "thin string" type (Kibble, 1976). The illustration 
will be extended to appropriately smooth perturbations of the metric. 

Roughly speaking, the method involves extending the Schwartz space of distributions to a much 
larger space, whose elements we will call generalised functions, within which the operations for computing 
the curvature can always be defined. A subspace of the generalised functions can then be defined whose 
elements correspond, in a many-one manner, to classical distributions. For those metrics whose curvature 
components lie in this subspace we can assign a well defined distributional curvature. This procedure 
circumvents the result of Schwarz quoted above by violating both the main conditions of Schwartz' 
theorem: in the space of generalised functions xS(x) ^ 0, and multiplication does not coincide with 
ordinary multiplication for continuous functions (although it does for C°° functions). These properties of 
generalised functions (which would usually be regarded as undesirable) are resolved when the generalised 
functions are mapped onto ordinary distributions, where this is possible. When this is done, xS(x) 
corresponds to the zero distribution; 
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2. An overview of Colombeau's generalised functions 



2.1. Smoothing distributions 

As we shall be concerned both with space-time, locally R 4 , and with a 2-dimensional plane transverse 
to the conical singularity, it will be convenient to work in R™ for general n. The basic tool will be the 
operation of smoothing distributions. Suppose <f> is a member of the space 2?(R") of test functions: 
smooth (i.e. C°°) C-valued functions on R" with compact support; and that 



or 



$(x)dx= 1. 
Given e > 0, we define 

so that $ e has a support scaled by e and an amplitude adjusted so that its integral is still unity. If 
/ : R™ — ► C is a function, not necessarily continuous, then by a smoothing of / we mean one of the 
convolutions 

/>) := y/(y + x)$(y)dy = f /(z)$(z - x) dz. 

/ e (x) := J /(y + x)$ £ (y)dy = J /(z)$ £ (z - x) dz. 

(Smoothed functions will always depend implicitly on $ as well as on the explicit arguments x and e.) 

Smoothing is defined in the same way for distributions, but with some notational changes. A distri- 
bution R is regarded as a C-valued functional 

V(R n ) 3 4> ^ (R, <f>) G C 

on the space T>(R n ) of test functions, and the convolution is defined by 

£(x) = (£,$(. -x)). 

i? e (x)-(i?,$ e (. -x)). 

An intuitively plausible procedure for defining the product Rf of, say, a distribution R and a dis- 
continuous function / would then be to define the action of the product on a test function ip by first 
defining the corresponding action of the product of the smoothed quantities R e and f e , and then taking 
the limit as the smoothing is made progressively finer, with e — > 0: 



{Rf,il>) = lim J 4(x)/ e (x)V(x)dx. 



For example, if R = 8, the Dirac <5-function, and / is the Heavisidc function, then this prescription yields 
the attractive solution 5f = 8/2, provided that $(— x) = 3>(x). If one considers products involving more 
complicated distributions, then the dependence of the answer on the nature of $ becomes more detailed. 
For example, writing x~ l for the distribution defined by taking the Cauchy principal value in integrals, 
x~ 1 S = kS' where A: is a ^-dependent constant. 

A further complication arises in the context of relativity, where one would like the results of multipli- 
cation to be invariant under C°° coordinate changes. The underlying problem here is that the operation 
of convolution depends on the linear structure of R 4 and so is not invariant under such changes. This 
suggests that smoothing by convolution should be regarded as a special case of a more general sort of 
smoothing, invariant under coordinate changes. 

Colombeau's definition of generalised functions therefore starts from a space £m(R") of functions 
depending on both the position x and a smoothing kernel <f>, on which are imposed conditions that 
reflect the special case of the smoothing of a distribution, but which are more general. Multiplication 
is defined on these pointwise. This produces a space with some of the required properties, but in which 
multiplication does not coincide with the ordinary multiplication even for C°° functions. This is rectified 
by defining an equivalence relation on £m(R") and passing to a space Q of equivalence classes. At the 
expense of minor complications, the whole of the construction may be applied to a domain which is a 
subset of R", in which case we denote the space of functions Em (O) and the space of equivalence classes 
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2.2. Specification of the algebra 

We now give a precise definition of the space of generalised functions Q(R n ). We shall use standard 
multi-index notation 

i = (ii, . . . ,i n ), |i| = ii H h i n - 

so 

X X\ . . . x n 

and 

D'= 3111 .. 
dxi 11 . . . dx n %n 

An essential role in defining both the basic space of functions £m(R") and the equivalence relation 
used to define the final generalised functions is played by a classification of smoothing kernels into subsets 
A q (R n ). 

Definition 1. For q e N we define .4 9 (R n ) to be the set of functions $ e V(R n ) such that 

(i) / RB $(x)dx=l 

(ii) J Rn 4>(x)x i dx = Vi e N™ such that |i| < q 
We then define £m(R™) in two steps. First, we set 

£(R n ) = { i? : A x R" -» C | x ^ R(<&, x) is C°° } 

This is an algebra (a vector space furnished with a multiplication) under the operation of pointwisc 
multiplication 

(i?S)($,x) := R($,x)S(&,x). 

The derivatives of these functions with respect to x will be denoted by D Q i?($, x). 
The second step is as follows: 

Definition 2. The subalgebra £m(R") (functions of moderate growth), of £(R n ) is defined to be the set 
of functions R such that for all compact K C R™ and for all a € N™, there is some N € N such that: If 
$ e An, 3c, 7] > such that 

\D a R($ £ , x)| < ce-^ (sei^, < e < rf) 

It is easily verified that, for any distribution T with support in R™, the corresponding smoothed 
function T is in £m(R™)- 

Finally, we need to define the equivalence relation that will give a multiplication that coincides with 
the usual one on C°° functions. The key idea is to note that there are two different ways in which ordinary 
functions / can be mapped into elements of £m(R"): wc can (for any continuous /) form by smoothing 
the element 

J2/(x,$) = y/(y + x)$( y )dy 
or (for C°° functions / only) we can form the element 

5 / (x,$)= / (x). 

Multiplication of the Sf trivially coincides with ordinary multiplication of the functions /, whereas this 
does not hold for the Rf. Colombeau therefore defines an equivalence relation that identifies Sf with 
Rf. This ensures that for C°° functions there is a single mapping of functions into generalised functions 
(equivalence classes) , which coincides with the mapping by smoothing used for distributions, and on which 
multiplication coincides with ordinary multiplication for C°° functions. The equivalence is defined by 
means of an ideal W(R") (i.e. a subalgebra such that if R e Af(R n ) and S e £ M (R") then RS e W(R")) 
so that elements differing by a member of Af(R n ) are identified. 
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Definition 3. The ideal Af(R n ), of Em (R n ) is defined to be the set of functions R such that for all 
compact K C R" and for all a £ N", there is some N £ N and some increasing and unbounded sequence 
{j n } such that: If $ £ *4 9 , for q > N, 3c, rj > such that 

\D a R(^,x)\<ce^- N (x £ < e < r?) 

For the smoothed functions used here, 7 g will be simply g; the more general form is used so as to 
ensure invariance under coordinate changes. 

The key result is then that, for C°° functions /, Rf — Sf is in Af(R n ). 

Definition 4 (Generalised functions). The algebra of generalised functions is defined to be 

£m(R") 



G(R n ) = 



W(R") 



The generalised function corresponding to the product of two distributions T and U , or of a distri- 
bution T and a continuous function /, is then taken to be the equivalence class of the product TU, or 
Tf respectively. 

2.3. Partial correspondence with distributions 

In order to associate a generalised function with a distribution we need to give a meaning to the expression 
/ G(x)0(x) dx with G a generalised function and (j) a test function. Colombeau therefore introduces a 
concept of the integral of a generalised function. 

Given a generalised function G = [i?( < &, x)], given as the equivalence class of an element R £ £m(R"), 
the quantity 



i?($,x) dx 



is a function of <& that depends on the choice of R. More specifically, it is a function p : A\ — > C satisfying 
the condition: 

there is some N £ N such that: If $ £ An, 3c, 77 > such that 

|p($ e )| < ce~ N (0 < e < rf) 

Definition 5. The algebra Em is defined to be the set of functions p satisfying the above condition. 

Different choices of the representative R for G result in functions p differing by a function in a 
subalgebra I playing the same role as Af(R n ): 

Definition 6. The ideal X, of Em is defined to be the set of functions R such that there is some N £ N 
and some increasing and unbounded sequence {7^} such that: If $ £ A q for q > N, 3c, 77 > such that 

|p($ 6 )| < cei"- N (0 < e < n) 

Thus if we take the quotient of Em by 1 we will obtain an equivalence class that is independent 
of the representative R, and that can therefore be regarded as a value of the integral of the generalised 
function G. Whereas the integral of an ordinary complex function is a complex number, an element of 
this quotient is called a generalised complex number. 

Definition 7 (Generalised Numbers). The algebra of generalised numbers is defined to be 

For each ordinary ordinary complex number z we can form the generalised complex number z defined 
as the equivalence class [p z ] of the constant function p z {<&) = z. It follows then that z is associated to z, 
in the following sense: 
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Definition 8. We say that z € C is associated to the classical number z E C (written as z h z) if there 
is some representative p € £m of z such that for $ e A q with g large enough, 

lim p($ e ) = z. 

Definition 9. We say that , z 2 € C are associated to each other if and only if z\ — z 2 h e C 

Note that zhz does not necessarily imply z — \p z \. It is, however, this weaker notion of association 
that turns out to be the important one in assigning corresponding distributions to generalised functions. 
To do this, we first make the following: 

Definition 10 (Weak equivalence). We say that G\, G 2 € Q(R n ) are associated to each other (written as 
G\ ~ G 2 ) if and only for each <PeP, 

/ (Gi(x)-G 2 (x))*(x)dxhOeC. 

Then we say that a generalised function G corresponds to a distribution T if G w T. For physical 
purposes, we are interested in those generalised functions, defined by replacing the terms in the defini- 
tion of the Riemann tensor components by smoothed functions and distributions, which correspond to 
distributions. 



3. Smoothing the cone 

In this section we will look at the curvature of the spacetime given by the metric (1). In this example the 
singular part of the curvature arises because of a conical singularity in the 2-dimensional plane transverse 
to the axis. We will therefore examine the curvature of the 2-dimcnsional metric 



ds 2 



dr 2 +A 2 r 2 d<f 2 . 



\A\<1. 



(2) 



Since this is the metric of a surface the above example has the advantage that the curvature may be 
given by a scalar quantity, the Gaussian curvature, rather than having to use the curvature tensor. This 
simplifies the notation and the details of the smoothing but plays no essential role in the calculation. We 
explore the implications of this further in the final section. 

We want to think of the cone metric as being defined on the whole of R 2 so that we start by 
expressing (2) in coordinates which are regular at the origin. In Cartesian coordinates (x = r cos <j), 
y = r sin <f>) the metric is given by 



g a b 



A 2 )S ab 



A 2 )h ab 



where 



h a b = 



y 



x A + y 2 

2xy 
x 2 + y 2 



2xy 



Note that h a b represents the 'singular' i.e. the non-smooth part of the metric. 

We now wish to regard the components of the metric as elements of Em (R 2 ) so that for $ e ^4i(R 2 ), 



9ab e (x,y) 



/ g a b(u,v)$ e (u 
JR 2 



X, v - 



y) ditdii, 



are smooth functions which we regard as the components of a smoothed metric. We can therefore calculate 
the Ricci curvature R e of this metric which may also be regarded as an element of £m(R 2 )- In this section 
we show that the curvature is given by a delta function in the sense that [-RVff] ~ 47r(l — A)5^ 2 \ 



Remark W. e shall regard the delta function as a scalar density, and the functions in D as scalars. Thus 
a factor of \Jg is inserted to make the left hand side a density. 
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It is important to note that although the calculation is carried out using the Cartesian components 
of the metric for convenience, the definitions of the Colombeau algebra ensure that the result does not 
depend upon the choice of coordinates used provided they are related by a smooth transformation. 

In order to calculate g^b e we nrs t note that constant functions are unchanged on smoothing so that 



g abfL {x,y) = i(l + A 2 )(U + i(l- A 2 )h abe 



(3) 



We next note that 



2 2 
X & - If 



2xy 



x 2 + y 2 x 2 + y 2 

So that we may calculate h a be by smoothing the complex valued function h(x, y) = e 21 * to obtain 



h € (x,y) = f h(u + x, v + y)$(u/e, v/e) dudv 
6 Jr. 2 



In order to calculate the integral we express it in terms of polar coordinates (r, <f>) and expand $ in terms 
of circular harmonics. Writing 



$(r, 



in<f> 



neZ 



we have 



where 



1 f 00 

~h e (r,4>) = -J2 In{r',r,<P)^n{r'/e)r'dr' 



In= [ 

Jo 



nez J ° 



(re 1 * + r'e 1 * )" 



d<f)'. 



r 2 + r' 2 + 2rr' cos(</> — <p') 

Setting w = e 1 * and z — c 11 * , we may integrate out 0' by considering the integral of the complex function 



T(z) 



r z + rw 



rz + r'w 



wz 



around the circular contour k : \z\ = 1. 

If r' < r, a single pole occurs at z = —r'w/r which has a residue of —\{—r'/r) n {\ — r' 2 / r)c^ n+2 ^ . 
If n — — 1, a single pole will occur at z — which has a residue of — i(r'/r)e I( ^. 
If n < -2, a |n|-fold pole will occur at z = which has a residue of i(-r'/r)"(l - / 2 /r)c l(n+2)l *. 
Using the residue theorem, we find that 



I n (r',r, 4>) 



So if we expand h e (r, <j>) as 



2tt [I - 
-2tt ( 1 



r.) e i(«+2)0 if / < r and n > 0, 
c i(«+2)0 if / > r an d n < -2, 



27r min{r'/r, r/r'je 1 * 
I 



if n = —1, 
otherwise. 



neZ 



we obtain 



27r(-l)"/ r/e (^) n (l-4^)$ n _ 2 (r / )r'dr', if n > 2, 

//„,(» = <( 27r^ /e (^)$„ 1 (r')r'dr' + 2^j^(^ 7 ) dr', ifn = l, 

^ -27r(-l)» /~ (^)" _2 (l - 4^) $n-2(r')r'dr', if n < 0, 

We next obtain estimates for the r dependence of H n ^(r). We will be interested in the behaviour both 
for small and large r/e. Let R = sup{r | |$ n (r)| > for some n € Z} then we consider exterior case 
r > ei?o first. 
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Case 1 (r > ei? )■ ■ 

In polar coordinates the condition that $ e A q may be expressed as 

/>oo 

/ $o(r)rdr = 1 
Jo 

r r ^ n (r)dr = 0, H< C <<Z 
J Q c + n is even 

Fixing r > and choosing e < r/i?o, we find that by (4) 

poo 2 /-''A 

i?2. e (r) = 1 - 2tt / r'$ (r') dr' - 2vr— / r' 3 $ (r') dr' 

Jr/e r JO 

Note that the first integral gives zero because r/e > R . So that 

H 2 ^r) = l + 0(e 2 /r 2 ) 
However if $ e *4 g with q>2, we may use (5) to express 77 2 , e as 

/■oo 2 /■oo 

ff 2 ,eW = 1 - 2tt / r'$ (V) dr' + 2tt— / r' 3 $ (r') dr' 

J r/e ^ J r/e 

in which both integrals vanish for e < r/R . A similar calculation shows that for n — 1 and for n > 2 

H n , e (r) = 0(e n /r n ) 

with a coefficient that vanishes if $ £ „4<j with q> n. 
On the other hand, by the definition of i?o 

ff„, £ (r) =0 for n < . 

Putting these results together we have in this case for $ e ^4 9 , 



(4) 
(5) 



h e (r,<P)=e*« + 0(^ 



Case 2 (r < eRq). In the following Ci, C2 etc. denote positive constants. 
For n > 2 or n < — 3 we use the fact that |<&| < K = sup |$| to obtain 



\H n J < 2ttK- 



1 



1 



n n + 2 



0(r7 



2 /,2\ 



The cases n = —2, — 1, 0, 1 are more delicate. Here we expand H n ^ e in a Taylor series (with remainder). 
For the case of H , e we have 



H 0te (r) = 2*J~ (l-^) *- 2 (r')dr' 
-2tt 



and so 



9 9ir 
$_ 2 (r')dr' + — [2$_ 2 (r/e) 



flo, e (0) = C\ £ 

^0,e(0) - 



*C(o) = § 



giving H 0t<i (r) = C3 + 0(r 2 /e 2 ). A similar method gives 

H_ 1>e (r)=0(r/e), 
H_ 2 , e (r)=0(r 2 e 2 ). 

Therefore in this case we have 

~h e (r, 0) = a + ai V -c^ + a 2 r -e^ + 0(r 2 /e 2 ) 

where a n are constants. 

The smoothed metric may now be expressed by (3), as 



h a b 



2 2 
x — if 



2xy 

x 2 +y 2 x 2 + y 2 
2xy y 2 — x 2 
x 2 + y 2 x 2 + y 2 



O 



(x 2 + y 2 )' 



for r > ei? , and 



hab 



fo + -(fox + foy) fo + -(fox + foy) 
e e 

02 + -{fox + foy) -fo - -(fox + foy) 



O 



x 2 + y 2 



for r < eR . 

We may now calculate the the Ricci scalar R € of the smoothed metric. 



Re 



0(1 /e 2 ) ifr<eRo 
C>(fS) ifr>ei? 



Note that if in Definition 3 we only consider compact sets K not containing the origin, then the second 
case in the above equation holds for small enough e, and R e then satisfies the conditions for membership 
of N. This can be interpreted as meaning that the curvature is concentrated at the origin. 

We now show that [Ry/g\ w 47r(l - A)S^ 2 \ It is sufficient to show that for each * e V(K 2 ), 

lim / R cy /f e (x,y)V(x,y)dxdy = 4tt(1 - A)*(0,0) 

for all <3? G A m for some m € N + , where if = supp \&. 

By the mean value theorem we may express the left hand side as ^(0, 0)2i + 12 where 



h 
h 



= / ReVthdzdy 

JK 

f _ Q^J 



for some £ € [0, 1]. Letting 



I'd* 



B e = {(x,y)GK 2 \ (x 2 +y 2 fl 2 KeRo) 



we may write 



U2I < c 



R t ^g e rdxdy 



+ c 


L 







R e ^/g e rdx dy 



For the first integral 



and for the second integral 



/ R e y/%rdxdy 

JB r 



< C 4 



(ei?o) 3 



where Rk is the maximum radius of K. 
This gives 



|/ 2 1 <C 6 e 



To calculate I\ we let D = { (a;, y) \ x 2 + y 2 < [i 2 } be a disc such that D C K and split up the 



integral into two parts. 



For the first integral we have 



h = / J R e y^Tdxdy+ / E eV /^dxdy 

JK-D JD 

[ R e y/g~ £ dxdy = 0(e«+ 1 /^+ 1 ) 

JK-D 



So that in the limit as e — > the only non-zero contribution comes from the integral over the disc. To 
calculate this we apply the Gauss-Bonnet theorem to convert it into an integral around the boundary. 



- / R e i/g~ e dxdy = 2ir - / Kg e ds 

Z JD JdD 



(6) 



where Kg e is the geodesic curvature of the smoothed metric. Since D lies in the exterior region for small 
enough e we get 



/ - 

JdD 



K 5f ds = 2nA + O 



9+1 



so that 



Therefore 



as claimed. 



I 1= 4tt- in A + O 



7FPT 



lim(*(0,0)/i +I 2 ) = 4tt(1-A) 

e—>0 



4. Perturbations 

If we introduce a perturbation on the metric that is small enough, we still obtain the same conclusion. 
Let us suppose we introduce an r 2 perturbation so that the metric takes the following form. 

ds 2 = dr 2 + A 2 r 2 (l + kr 2 )d<j) 2 A,k = constant 

In Cartesian coordinates, the metric may be written as 

g ab = + A 2 (l + kr 2 ))6 ab + 1(1 - A 2 )h ab - hzA 2 l ab 



with the additional perturbation term being 



lab — 



x — y 2xy 
2xy y 2 — x 2 



Since a smooth function may be regarded as a generalised function with no $ dependence, an element of 
£(R 2 ) equivalent to the smoothed metric may be defined by 

g^ bt {x,y) = ^{l + A 2 +A 2 kr 2 )S ab + ^(l-A 2 )K b - ^A 2 kl ab 



We now consider the Ricci scalar of this metric. The main difference between this case and a flat 
cone is that the curvature does not vanish over away from the origin. In fact 



R 



-2k 



3 + 2kr 2 
(1 + fcr 2 ) 2 ' 



Using the smoothed metric, 



Re = 



0(l/e 2 ) if r < ei? 

R + O (fS) if r > eR 



We might expect to find that, in the sense of generalised functions, [-RVff] ~ R\/9 + 4tt(1 — A)5^ 2 \ 
We can verify this by performing a calculation similar to the one above for the flat cone, but by replacing 
Ry/g by Ry/g — R\/g throughout. 

The calculation of h is unchanged. To calculate h we write the integral as 



h 



ReVfedxdy - / Ri/gdxdy+ / {R e xjg\ - R^/g) dx dy 

D JD JK-D 



A l + 2k^ 2 

= 4tt - At: A- , 9 , 1/9 + 47rA 

(1 + ky 2 ) 1 / 2 



1 + 2fc^ 2 



(l + fc M 2 ) 1/2 



O 



9+1 



where we have applied the Gauss-Bonnet theorem (6) to evaluate the first integral, and explicitly com- 
puted the second. Therefore 

lim(*(0,0)/i +I 2 ) = 4tt(1 - A). 



and hence 



[Ry/$\**Ry/9 + *K<l-A)5<> 



2) 



It is clear from the above that replacing kr 2 by a more general perturbation of the form k{r)r 2 where 
k(r) is a smooth function would give the same result. Indeed any perturbation for which l a b is a smooth 
function would give a similar result. 



5. Conclusion 

We have shown that the curvature of a 2-dimensional cone at r = 0, when regarded as a generalised 
function is weakly equivalent to the delta distribution and the conclusion remains the same for a suitably 
perturbed metric. To simplify the presentation we have confined our attention to the two dimensional 
case however similar results hold for the case of the the 4-dimensional cone-like space-time (1). In the 
four dimensional case one must use a more general smoothing function with both t and z dependence 
to calculate the transverse components of g^b e - However after integrating out in the t and z directions 
the calculation is much the same. One may then calculate the curvature tensor R a bcd of the smoothed 
metric and regard the components as elements of Sm (R- 4 )- The only terms that might correspond to 
distributions are those which are the sectional curvatures of the 2-surface transverse to the string. By 
applying the methods of §3 but using a holonomy argument rather than the Gauss-Bonnet theorem to 
calculate the final integral one may show that these terms regarded as generalised functions are weakly 
equivalent to delta functions, thus confirming the results of Vickers (1987). 

Currently work is in progress to extend these calculations to more general four dimensional space- 
times in general relativity which admit quasi-regular singularities. One result would be to derive a new 
definition of the mass of a string per unit length, which Geroch and Traschen showed was ambiguous 
when defined by limiting sequences of metrics. 
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